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ABSTRACT 

This paper supplements Jiang et al.' ('2003'), who studied 172 M31 globular clusters (GCs) and globular cluster 
candidates from Battistini et al. (1987) on the basis of integrated photometric measurements in the Beijing- 
Arizona-Taiwan-Connecticut (BATC) photometric system. Here, we present multicolor photometric CCD data 
(in the BATC system) for the remaining 39 M31 GCs and candidates. In addition, the ages of 35 GCs are 
constrained by comparing our accurate photometry with updated theoretical stellar synthesis models. We use 
photometric measurements from GALEX in the far- and near-ultraviolet and 2MASS infrared JHKs data, in 
combination with optical photometry. Except for two clusters, the ages of the other sample GCs are all older 
than 1 Gyr. Their age distribution shows that most sample clusters are younger than 6 Gyr, with a peak at ~ 3 
Gyr, although the 'usual' complement of well-known old GCs (i.e., GCs of similar age as the majority of the 
Galactic GCs) is present as well. 

Subject headings: galaxies: individual (M31) - galaxies: star clusters - galaxies: stellar content 



1. INTRODUCTION 

The process of galaxy formation and evolution ranks among 
the most importan t outstanding problems in astrophysics (e.g., 
iPerrett et alj2002 !). One way to better understand the underly- 
ing questions is by studying globular clusters (GCs). GCs are 
often considered the fossils of the galaxy formation process, 
since they tend to form i n the very ear ly stage of their host 
galaxy's evolution (Bar mby et alJ l2000). A GC is a densely 
packed, gravitationally bound, roughly spherical system of 
several thousand to about one million stars. They can be ob- 
served out to great distances, which implies that they can be 
used to study and probe the properties of extragalactic sys- 
tems. The most distant GC systems studied to date are located 
in the Coma Cluster; their study has been facilitated by Hub- 
ble Space Telescope (HST ) Wide Field and Planetary Camera- 
2 (WFPC2) observations (Baum et al."1995'; 'Kavelaa rs et"an 
HbOO; Harris et al. 2000; Woodworth & Harris 2000). 

M31 (NGC 224), the Andromeda galaxy, is an early- 
type spira l galaxy (type Sb), located at a distan ce of ~ 
780 kpc (IStanek & GarnavichI IT998t iMacril IIOOTI) . It is 
the nearest and largest spiral galaxy in the Local Group 
of galaxies and has been th e subject of many GC studies 
and surveys. Hubbl^ (Il932h first discovered 140 GC can- 
didates characterized by Wpg < 18 mag. Subsequently , 
a number of studies dSev fert & Nassau 1 1 9451: iHiltned 1 1 9581: 
iMavall & EggeiJ 119531: IXron & Maval]ffT960 ') identified 



160 GC candidates in M31. IVetesnikI (11962) compiled 
the first major M31 GC catalog, containing about 300 
GC candidates and including UBV photometric data. The 
most extensive GC surveys have since be en published by 
ISargen t et ar('1 977h.lCram pton et al.l(| ' 
group (Battisti ni et al.l 11980, 198^ 



ICrampton et all (198 5*) and B attisti ni et allTl 980l 1 1 987ll 1 993h 



1985), and the Bologna 
19931). In particular . 



provided photometric data in either UBV or UBVR. These 
surveys were mostly based on visual searches of photo- 
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graphic plates and are fairly complete dow n to V - 18 mag 
(My ~ -6.5 mag) dFusi Pecci et al.|[l993h . although a num- 
ber of recent studies searched for fainter GCs in M31 (e.g., 
Mochejska et al. 1998; Barmby & Huchra 2001; Ki m et alj 
120071) . However. IGalleti et al.l (1200 6) showed that a signifi- 
cant number of class-D and E GCs with V > 17 are still to be 
confirmed (and hence the GC luminosity function is incom- 
plete), and that large surveys are needed to reach a complete 
sample of M31 GCs. Following on from the first extensive 
spectroscopic survey of M31 GCs by van den Bergh ( 1969), a 
signific ant number of authors (e.g., Huchra et al. 1982, 19911; 
i Dubath & Grillmair 1997; Federici et al. 1993; Jablonk a etalJ 
Il998'; Barmby et al. 200Q: iPeiTett et al.. ,2002;: ,Galleti et al.1 
2006; Lee et al. 20081 and references therein) embarked on 
studies of their spatial, kinematic, and metal-abundance prop- 
erties. The first comprehensive catalog including photomet- 
ric and spectroscopi c data for M31 GCs was assembled by 
iBarmbv et al](l2000h . The Revised Bologna Catalog (RBC) of 
M31 GCs was recently published by Galleti et al. (200%and 
has since been revised a number of times (Galleti et al. 20051 
2006, 2007). In the primary catalog (Galleti et al. 2004), all 
known M3 1 GCs and candidates were compiled based on a lit- 
erature survey, leading to a total of 1 164 entries including 337 
confirmed GCs, 688 GC candidates, and 10 o bjec t s with unde- 
termined classification. In addition. Galleti et al.l (12004) iden- 
tified 693 known and candidate GCs in M3 1 using the 2MASS 
database and included their 2MASS JHKs magnitudes. The 
latest RBC (V3.5) was updated on March 27, 2008, and in- 
cludes t he newly discove red st ar clusters fromlMackev et"al] 
(2006»), iKim et al.l (l2007l) . and iHuxor et al.1 (2008^ In total, 
1567 GCs and GC candidates (509 confirmed GCs and 1058 
GC candidates) are known in M31; 421 former GC candi- 
dates turned out to be stars, asterisms, galaxies, Hii regions, 
or extended clusters. In addition, the RBC V3.5 includes pho- 
tometric data of M3 1 GCs and GC candidates in the far- and 
near-ultraviolet (FUV and NUV) from the Nearby Galaxies 
Survey (NGS) o f the Galaxy Evol ution Explorer (GAL EX) 
dRev et al.ll20"07h . Very recently, ,Caldwell et all ( l2009t) pre- 
sented a new catalog of 670 Ukely star clusters in the field 
of M3 1 , all with updated high-quality coordinates accurate to 
0.2", based on images from either the Local Group Survey 
(lMasseyll2006h or the Digital Sky Survey. 
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An accurate and reliable analysis of integrated stellar pop- 
ulations (such as star clusters) is key to understanding the 
formation and evolutionary process in galaxies. By means 
of comparisons of integrated populations with models of ho- 
mogeneous stellar systems, i.e., simple stellar populations 
(SSPs), recent studies have achieved some success in de- 
termi ning ages and mas ses for extragalactic star clusters 
(e.g..lde Gri is et al.'2003a'W.l Bik et al.l200llMa et al.l2006l: 
iFan et alJl20 06: Ma et al. 20071 

iMa et all (1200 6) and Fan et alj (2006') derive d age esti- 
mates forM31 GCs by fitting SSP models ( Bruzual & CharlotI 
I2003L henceforth BC03) to their photometric measurements 
in a large number of intermediate- and broad-band passbands 
from the optical to the near-infrared (NIR). For instance, 
iMa et al.1 (l2b07i) constrained the age of the M31 GC S312 
(B379), using multicolor photometry from the NUV to the 
NIR, to 9.5+^ '^ Gyr. S3 12 (B379) is among the first extra- 
galactic GCs for which the age was estimated accurately us- 



:iy 

ing main-sequence photometry, i.e., iBrown et al.l (120041) es- 
timated its age at lO^j'^ Gyr. This was based on their anal- 
ysis of the cluster's color-magnitude diagram (CMD) below 
the main-sequence turn-off" (MSTO) using extremely deep im- 
ages obtained with the Advanced Camera for Surveys (ACS) 
on board the HST. They performed a quantitative compari- 
son of their resolved stellar photometry with the isochrones 
of VandenBerg et al. (2006). 

In this paper we first describe our new observations and the 
relevant data-processing steps, as well as the complementary 
data used from the literature (®. In ^we quantitatively 
compare the spectral energy distributions (SEDs) of the GCs 
in our sample with the galev SSP models. Finally, our results 
and a summary are presented in ^ 

2. DATABASE 

2.1. The sample 

The GC sample u sed in this pa p er wa s taken from the 
Bologna catalog of Batt istini et al.l ( 119871) . which contains 
827 M31 GCs and GC can didates. In addit ion, our sam- 
ple also supplements that oflJiang et al.l (2003b. who studied 
172 GC candidates from 'Battisti ni et all (Il987l) on the ba- 
sis of integrated photometric measurements in the Beijing- 
Arizona- Taiwan-C onnecticut (BATC) photometric system. In 
iJiang et all (l2003l). all GC c andid ates of classes A and B (353 
objects) in iBattistini et all (Il987h (i.e., their Table IV) were 
adopted as their original sample. However, of these only 223 
objects are in their observed CCD fields. They also noted 
that B007 is a gala xy, and B05 5, B132 and B147 are vk- 
tually stars (.Barmby et al.l l2000h. These four objects were 
therefore not included in ljiang et al.l (l2003h 's final sample. In 
summary, 219 class-A or B GCs were observed bv lJiang et al.l 
(|2003), of which 47 were excluded because of missing pho- 
tometric measurements in some filters. In this paper we ana- 
lyze these 47 GC candidates on the basis of newly observed 
data in the BATC photometric system combined with GALEX 
FUV/NUV photometry, broad-band UBVRI, and NIR JHKs 
(2MASS) data. However, we did not manage to obtain accu- 
rate photometric measurements for a number of objects be- 
cause of either the dominance of a nearby very bright ob- 
ject (B095, B176, and B202), the GC candidate being very 
faint and superimposed on a very high background (Bl 19 and 
B324), or the GC candidate being located near M32 (B124) or 
NGC 205 (B331), both also resulting in a very high contribu- 
tion. In addition, object B330 is faint and located very close to 
a brighter object, rendering accurate photometry impossible. 
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Fig. 1 . — Spatial distribution of the GC candidates in M3 1 . Circles and plu s 
signs represent the samples discussed in this paper and bv lJiang et alj 120031) . 
respect i vely. T he large ellipse is the M3 1 disk/halo boundary as defined by 
IRacind 09911) : the two small ellipses are the D25 isophotes of NGC 205 
(northwest) and M32 (southeast). 



Thus, here we analyze the multicolor photometric properties 
of 39 GC candidates. Figure 1 shows the s patial distributions 
of both our sample GCs (circles) and the I Jiang et aP (l2003l) 
GCs (plus signs) across the M3 1 field. 

2.2. Observations and data reduction 

To obtain photometric measurements in the B ATC photo- 
metric system for the 39 GC candidates for which J iang et al] 
(l2003l) did not obtain photometry in a number of filters, we 
re-observed the objects. The BATC photometric system uses 
a Ford Aerospace 2048x2048 CCD camera with a pixel size 
of 15 fj.m, mounted at the focus of the 0.6/0. 9m f/3 Schmidt 
telescope at Xinglong Station (National Astronomical Obser- 
vatories of China; NAOC). The CCD field of view is 58'x58', 
with a pixel size of 1.7". The multicolor BATC filter system 
includes 15 intermediate -band filters covering the wavelength 
range from 3300A to 1 fim. These filters were specifically 
designed to avoid contamination from the brightest and most 
variable night-sky emission lines. The CCD camera is not 
sensitive at the shortest wavel engths covered by the BATC fil- 
ters. For this reason, neither IJiang et alj (l2003h nor we used 
the two bluest filters (a and b) for our observations. Finding 
charts of the sample GCs and GC candidates in the BATC g 
band (centered at 5795A), obtained with the NAOC 60/90cm 
Schmidt telescope, are shown in Fig. 2. 

Thirteen hours of imaging of the M31 field of lJiang et alj 
(2003) were obtained through the usable set of 13 
intermediate-band filters from November 15, 2003 to Decem- 
ber 13, 2003. Bias subtraction and flat fielding using dome 
flats were done with the BATC automatic data-reduction soft- 
ware, PIPELINE I, ori ginally developed for the BATC Multi- 
color Sky Survey ( Fan et alj |1996': 'Zheng et al. 1999). The 
dome flat-field images were taken using a diff'user plate in 
front of the Schmidt telescope's corrector plate. This flat- 
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fielding technique was verified using photometry obtained for 
other galaxies and spectrophotometric observations (see, e.g., 
Fanet al."1996"; 'Zheng et al." 19991; IWu et al.ll2002t lYan et al.l 
2000; Zhouetal. 2001, 2004). Spectrophotometric calibra- 
tion of the M31 images was done by observations of four 
F-type subdwarfs, HP 19445, HP 8 4937, BP +26°2606, 
and BP +17°4708 (lOke & Gunnlll983l) . Our magnitudes are 
therefore defined in the spectrophotometric AB magnitude 
system (i.e., the Oke & Gunn monochromatic system), 

mBATC = -2.51ogFv- 48.60, (1) 

where Fy is the appropriately averaged monochromatic flux 
in units of erg s"' cm"^ Hz"' at the effective wavelength of 
the specific pas sband. In the BATC system Fy is defined as 
(lYan et al.ll2000h 

/d(logv)/,r, 

Fv = — p , (2) 

J d{\ogv)ry 

which relates the magnitude to th e number of photons de- 
tected rather than to the input flux (iFukugita et al.|[T99 ^. In 
Eq. (2), Fy is the system's response and fy the object's SEP. 
Spectrophotometric calibration of the M3 1 images using the 
Oke-Gunn sta ndard stars was done during photo metric nights 
(see for details lYan et al.l2000l;IZhou et al.ll2001h . Using these 
standard-star images, we iteratively obtained atmospheric ex- 
tinction curves and the variation of these e xtinction coeffi- 
cients as a funct ion of the time of night (cf. lYan et al.ll2000l; 
IZhou et al.ll200Th . 

mBATC = mi„,t + {K + ^K{\5T)]X + C, (3) 

where X is the airmass and [K + AK(UT)] the time-dependent 
extinction term. The instrumental magnitudes (;«inst) of se- 
lected bright, isolated, and unsaturated stars on the M3 1 im- 
ages observed on photometric nights can be readily trans- 
formed to the BATC system (»jbatc)- The calibrated magni- 
tudes of these stars were then used as secondary standards to 
uniformly combine images from calibrated nights with their 
counterparts observed on non-photometric nights. Table 1 
lists the parameters of the BATC filters and the observational 
statistics; column 6 provides the scatter, in magnitudes, for the 
photometric observations of the four primary standard stars in 
each filter. 

2.3. Integrated photometry 

For each M31 GC cand idate we used the phot routine in 
PAOPHOT (IStetsonlll987i) to obtain the integrated photome- 
try. To avoid contamination from nearby objects, we adopted 
an aperture of 10.2" diameter, corresponding to 6 pixels. In- 
ner and outer radii for background determination were taken 
at 8 and 1 3 pixels from the GC center. Given the small aper- 
ture used for the GC observations, aperture corrections were 
determined as follows. We used isolated stars to determine 
the magnitude difference between diameter of 6 pixels and 
the fully integrated magnitude of these stars in each of the 1 3 
BATC filters used. The SEPs for our sample of 39 GCs and 
GC candidates were then corrected for the filter-specific dif- 
ferences, and these values are given in Table 2. Columns 2-14 
give the magnitudes in the 13 BATC passbands used for our 
observations. For each object the second line lists the Icr un- 
certainties in magnitude for the corresponding passband. The 
errors for each filter are given by PAOPHOT. The magnitudes 
of B129 in the c and d filters, and that of B195 in p filter could 



not be obtained owing to low signal-to-noise ratios in these 
filters. 

2.4. GALEX, broad-band, and 2MASS photometry 

To estimate the ages of the M31 sample GCs and GC 
candidates accurately, we use as many photometric data 
poi nts covering a s wide a wavelength rang e as possible 
(cf. 'deGriis et al."2003bt lAnders etal.ll2004 . In addition, 
Kaviraj et al. (2007) showed that the combination of FUV 
and NUV photometry with optical observations in the stan- 
dard broad bands enab les one to effic iently break the age- 
metallicity degeneracy. IWorthevI (1 19941) showed that the age- 
metallicity degeneracy associated with optical b road-band 
colors is A age/AZ - 3/2 (see also iMacArfliure t al. 2004). 
However, de Jond (11996) showed that this degeneracy can 
be partially broken by adding NIR phot ometry to optica l 
colors, which was recently confirmed by IWu et al.l (l2005l) . 
ICardiel et al.l (l2003h found that the inclusion of an infrared 
(IR) passband can improve the predictive power of the stel- 
lar population diagnostics by ~ 30 times compared to using 
optical photometry alone. Since NIR photometry is less sen- 
sitive t o interstellar extinc t ion th an th e classical optical pass- 
bands, iKissler-Patig et al.1 (l2002h and lPuzia et"an (12002 ) also 
suggested that it provides useful complementary information 
that can help to disentangl e the age-metallicity degeneracy 
(a lso seelGalleti et al.| |2004l). 

iRev et al.1 (l200l published GALEX NUV and FUV pho- 
tometric data for 485 and 273 M31 GCs, respectively. The 
photometric data for 28 (NUV) and 17 (FUV) of our M31 
sample GCs in common is listed in Table 3. Again, for each 
object the second line lists the photometric uncertainties for 
the corresponding passband. The GALEX photometric sys- 
tem is calibrated to match the spectrophotometric AB system. 

To date, the study of M3 1 GCs has been largely based on th e 
excellent Bologna catalog (Bat tistini etakll 1 98 0'. 1 987, ' 1 9931). 
Updates to the original RBC were provided by Galleti et al] 
(2004) who take as t heir photometric reference the dataset of 
,Barmbv et alj(l2000l) in order to obtain the most homogeneous 
set of photometric measurements available. iBarmbv et alj 
(2000^) published optical and IR photometric data for 285 
M31 GCs (see their Table 3), obtained with the 4-Shooter 
CCP mosaic camera and the SAO IR imager on the 1.2m 
telescope at the Fred Lawrence Whipple Observatory. Photo- 
metric measurement s in the UBVRI bands were published by 
iBarmby et a l.' (2000") for most of our sample objects. There- 
fore, we pref erentially adopt the UBVRI measurements of 
Barm bv et al.l (120 00). For the remaining GCs we follow 
iGalleti et al.l (1200 4). who updated the Bologna catalog with 
homogenised optical (UBVRI) photometry collec ted from the 
most recent photometric references available. iGalleti et alj 
(2004) did not include the photometric uncertainties. Al- 
though we refer to the original papers, the uncertainties as- 
sociated with the same object but based on the use of differ- 
ent photome t ric sys tems are often very different. In addition, 
IGalleti et all ( 1200 4*) transform ed their UBVRI photometry to 
the reference system of Barmb vet al.l 1200 0) by applying off- 
sets derived from objects in com mon be tween the relevant 
catalog and the data set of .Barmby et alJ (I20b0i) . The mea- 
surements are therefore internally consistent, and referenc- 
ing the original uncertainties may be irrelevant. Therefore, 
we only adopted p hotometric uncertainties as suggested by 
iGaUeti et alJ (l2004l) . i.e., 0.05 mag in BVRI and 0.08 mag in 
U. In fact, these photometric uncertai nties do n o t affec t our 
results significantly, as we showed in iFan et al.] (l2006l) (see 




Fig. 2. — Finding charts of the sample GCs and GC candidates in the BATC g band, obtained with the NAOC 60/90cm Schmidt telescope. The field of view of 
each image is 11' x 11'. 
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th eir §4.3 for d e tails). 

iGalleti et aH (l2004l) identified 693 known and candidate 
GCs in M3 1 using the 2MASS database, and determined their 
2MASS JHKs photometr ic magnitudes (transfo rmed to the 
CIT photometric system) (lElias et al.l ll982. 1983). However, 
we need the original 2MASS JHKs magnitudes for our sam- 
ple GCs to compare our observational SEDs with the SSP 
models, so we rev ersed this tra nsformation using the same 
procedures. Since IGalleti et al.l (|2004) did not provide the 
2MASS JHKs uncertainties, we obtained the se by comparing 
the ph otometric magnitudes with Fig. 2 of I Carpenter et al.l 
(1200 ll) . They show the observed photometric rms uncertain- 
ties as a function of magnitude for stars brighter than their ob- 
servational completeness limits. We include the broad-band 
and 2MASS photometry (and the associated uncertainties) of 
the sample GCs in Table 3. We also list the new classification 
flags, following RBC V3.5 notation. From Table 3 we learn 
that B052 and B062 are classified as galaxies based on their 
radial velocities. We will therefore not estimate their ages be- 
low. 

2.5. Comparison with previously published photometry 

The BATC intermediate-band system can easily be trans- 
formed to the f/BVT?/ broad-band system. Zhou et al. (2003]) 
derived the relationships between these two sys t ems using 
standard stars from the catalo gs of lLandollI (Il983l 1 1992) and 
iGaladf-Ennquez et al.l (l2000h : 



niB = nid + 0.2201(m, - m,) + 0.1278 + 0.076, (4) 



mv = nig + 0.3292(m/ - m/,) + 0.0476 + 0.027. (5) 

To check our photometry we derived the magnitudes in B 
and V based on Eqs. (4) and (5). We transformed the mag- 
nitudes of our 39 GCs and GC candidates in the BATC c, d, 
e bands to B-band photometry, and BATC /, g, and /z-band 
measurements into V-band data. Fig. 3 shows a compari- 
son of our V and (B - V) photometry with pre viously pub- 
lished measurements of lBarmbv et al.l (120001) and lGalleti et al.l 
( l2004l) . The mean V magnitude and (B - V) color differ- 
en ces (in the sense o f this paper minus Barmbv et al] (l2000h 
or lGaUeti et aH ("2004)) are (AV) = -0.066 + 0.013 mag and 
(A(B - V)) = -0.040 + 0.017 mag, respectively, thus showing 
excellent agreement. 

2.6. Metallicities and reddening values 

To estimate the ages of our sample GCs accurately we 
required that our GCs have both independently determined 
metallicities and reddening values. We used t hree homoge- 
neous sources for spectroscopic metallicities dHuchra et al.l 
[1991; Barmby et al. 2000; Perrett et al. 2002*) and one refer- 
ence (Fan et al. 2008). 

iHuchra et al.i (il991i) obtained spectroscopy of 150 M31 
GCs and candidates with the Multiple Mirror Telescope 
(MMT). The system they used has a resolution of 8-9A and 
enhanced blue sensitivity. To obtain many of the strongest 
and most metallicity-sensitive spectral features of interest in 
the ultraviolet, they extended their observations to the at- 
mospheric cut-off at 3200A (see details in iBrodie & Huchra I 
[l990). The metallicities of these 150 objects were determined 
using six absorption-line ind ices from integrated cluste r spec- 
tra employing the method of iBrodie & Huchral (Il990h . 
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Fig. 3. — Comparison of our newly obtained cluster photom etry with 
previo us measurements by Bannby et al. (2000) (triangles) and Ga lleti et all 
(20041) (crosses). The dashed lines enclose ±0.2 mag in V and ±0.3 mag in 
B-V. 



iBarmbv eta"n(l2000h observed 61 M31 GCs and candidates 
spectroscopically using the Keck Low Resolution Imaging 
Spectrometer (LRIS) and the MMT Blue Channel spectro- 
graph. With Keck LRIS, they used a 600 ( mm ' grating with 
a 1.2A pixeL' dispersion from 3670-6200A, and a resolu- 
tion of 4-5 A. With the MMT Blue Channel, they used a 300 i 
mm"' grating with a 3. 2 A pixeL' dispersion from 3400-7200 
A, and a resolution of 9-1 1 A. They obtained the clu ster metal- 
licities on the basis of the IBrodie & Huchral(ll990h method as 
well. 

iPerrett et al.l (12002b determined metallicities for more than 
200 M3 1 GCs and candidates using the Wide Field Fibre Op- 
tic Spectrograph (WYFFOS) at the 4.2m William Herschel 
Telescope. Their spectral range covers ~ 3700-5600A using 
two gratings, one of which (H2400B, 2400 ( mm"') yields a 
dispersion of 0.8A pixel"' and a spectral resolution of 2.5A 
over the range 3700^500A, and the other (R1200R, 1200 
{ mm"') yields a dispersion of 1.5A pixel"' and a spectral 
resolu tion of 5.1 A over the range 4400-5600A. .Perrett et al.l 
(120021) cal culated 12 absorption-Une indices, again using the 
method of ^ Brodie & H uchra (1990). Through a comparison 
of the line indices with published M3 1 GC [Fe/H] value s fro m 
previous studies dBonoli et al.ll 19871; IBrodie & Huchra1ll990l: 
iBarmbv et al.l2000l) . thev found that the line indices of the CH 
(G band), Mg^i, a nd Fe53 lines best re presented their observed 
GCs. Therefore, iPerrett et al.l (120021) determined the metal- 
licities of their sample targets using an unweighted mean of 
these three [Fe/H] values. 

Using meta llicities fro m the hterature (iHuchra et al.lll99l1; 
IBarmbv et all 120 00; Perr ett et al.1 12002|) combined wifli flie 
RBC. iFan et al.l " il2008.) determined 443 reddening values and 
intrinsic colors, as well as 209 metallicities for individual 
clusters without spectroscopic observations. 

To use all metallicities as coherently as possib le we ranked 
the sou rces of M31 GCs metallicities, choosing IPerrett et al.l 
(I2OO2I) metallicities whenever available because of the large 
number of metallicity d eterminations . Me tallicities from 
Barmbv et al. (2000) and Huchra et al.l (Il991l) were preferred 
if Perrett et al. ( 2002) determinations were not av ailable. If 
spectroscopic metallicities were missing, we used iFan et al.l 
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( 120081) . Metallicities were not available for B089 and B226. 
As a consequence, we do not attempt to determine their ages 
(see details in §4). The final set of metallicities for the sample 
clusters is included in Table 4. 

For the r eddening value s of the s a mple G Cs we refer to 
iBarmbv et al. (.2000.) and iFan et alJ (l2008h . iBarmby et all 
(I2000h determined the reddening for each cluster using corre- 
lations between optical and IR colors and metallicity, and by 
defining various 'reddening-free' parame ters using the i r large 
database of multicolor photometry. Barmbv e t al.l ( l2000h 
found that the M31 and Galactic GC extinction laws, and 
the M3 1 and Galactic GC color-metallicity relations are sim- 
ilar They estimated the reddening to M31 objects with spec- 
troscopic data using the relationship between intrinsic opti- 
cal color and metallicity for Galactic clusters. For objects 
without spectroscopic data they used the relationships be- 
tween the reddening-free parameters and certain intrinsic col- 
ors based on Galact ic GC data . Follo wing the methods in 
IBarmby et aP (l2000h . iFan et al.1 (l2008h (re-)determined red- 
dening valuesfor 443 clusters and cluster candidates. We 
choose lFan et alj ( l2008l) reddening values whenever available 
because their reddening values co mprise a homogeneou s data 
set and they are larger than those of Barmby et al.l (12000) . The 
reddening values for the sample clusters are listed in Table 
4. The values of extinction coefficient are obtaine d by in- 
terpola ting the interstellar extinction curve of iCardelh et al.l 
([T989h . 

3. AGE DETERMINATION 
3.1. Stellar populations and synthetic photometry 

The most direct method to constrain the ages of different 
stellar populations involves comparing the observed luminos- 
ity levels of the MSTOs. Unfortunately, this approach is lim- 
ited to the nearest GCs, where individual stars can be re- 
solved and measured down to a few magnitudes fainter than 
the MSTO. Even in M31, the nearest large spiral galaxy, 
the MSTO is only reached for one GC (S3 12) ( also see 
iBrownet al.1 120041: iRev et al.1 1200^ iMa et al.lll)07h. H ow- 
ever, since the pio neering work of iTinslevI (1 196811 1972 1) and 
ISearle et"an (Il973h evolutionary population synthesis model- 
ing has become a powerful tool for the interpretation of inte- 
grated spectrophotometric observations of gal axies as well as 
their components (see e.g. lAnders et al.ll2004l) . 

In evolutionary synthesis models, SSPs are modeled on 
the basis of a collection of evolutionary tracks of stars 
of different initial masses and a set of stellar spectra at 
different evolutionary stages. To estimate the ages of 
our sample GCs w e compare their SEDs with the galev 
SSP models (e.g., iKurth et a ll IT999t ISchulz etal.1 l200l 
lAnders & Fritze-vT~A lvenslebenl 2003h . The galev SSPs are 
based on the Padova isoch rones (w ith the most recent versions 
using the updated B ertelli et al.l {l994) isochrones, which 
include the thermal ly-pulsing asym ptotic giant-branch [TP- 
AGB] phase), and a Salpeterld 19551) stellar initial mass func- 
tion with a lower-mass limit of 0.10 Mq and the upper-mass 
limit between 50 and 70 Mq, depending on metallicity. The 
full set of models spans the wavelength range from 91 A to 
160 fim. These models cover ages from 4 X 10^ to 1.6 x lO"' 
yr, with an age resolution of 4 Myr for ages up to 2.35 Gyr, 
and 20 Myr for greater ages. 

Since our observational data consists of integrated lumi- 
nosities through a given set of filters, we convolved the theo- 
retical SSP SEDs with the GALEX FUV and NUV, broad- 



band UBVRI, BATC, and 2MASS JHK, filter response 
curves to obtain synthetic ultraviolet, optical, and NIR pho- 
tometry for comparison. The synthetic magnitude in the AB 
magnitude system for the /* filter can be computed as 



nii = -2.5 log 



-48.60, 



(6) 



where is the theoretical SED and (pi the response curve of 
the /* filter of the GALEX FUV/NUV, UBVRI, BATC, and 
2MASS JHKs photometric systems. Here, changes as a 
function of age and metallicity. 

3.2. Fitting results 

We use ax^ minimization test to examine which galev SSP 
models are most compatible with the observed SEDs, follow- 
ing 



mod "1 2 



(7) 



where OT™''(f) is the integrated magnitude in the /* filter of 
a theoretical SSP at age f, m™" represents the intrinsic inte- 
grated magnitude in the same filter, and 



2,2 
''"obs.i + ''"mod,/- 



(8) 



Here, cr\ . is the observational uncertainty, and cr^ . . is the 

' obs,? ' mod,! 

uncertainty a s sociat ed with the model itself, for the /* filter 
ICharlot et al.l d 19961) estimated the uncertainty associated with 
the term cr^^j , by comparing the colors obtained from differ- 
ent s tellar e volutionar y track s and spectral libraries. Follow- 
ing |WueLay (|2005|), |Ma.eLaL| ((20^ and lFanetaTJ(l2006h 
we adopt cr^ , . = 0.05. In fact, the values <t^ . . adopted do 

^ mod,? mod,; ^ 

not change the best fits, but only affect the;i'^ values. 

The GALEV SSP models include five initial metallicities, 
Z = 0.0004,0.004,0.008,0.02 (solai- metallicity), and 0.05. 
Spectra for other metallicities can be obtained by linear inter- 
polation of the appropriate spectra for any of these metallici- 
ties. In addition, if the metallicity of a cluster is poorer than 
Z = 0.0004, we only use the model of Z = 0.0004. The best 
fits to the SEDs of our GCs are presented in Fig. 4. 

4. RESULTS AND SUMMARY 

In the previous Section we determined the ages of 35 M3 1 
GCs and GC candidates. The results are listed in Table 5. The 
metalHcity of B089 and B226, and the reddening of B089 had 
not been determined previously. From Fig. 5, which shows 
the age distribution of the sample clusters (see also Table 5) 
we conclude that, except for two clusters, the ages of the other 
sample GCs are all older than 1 Gyr Most sample GCs are 
younger than 6 Gyr, with a peak at ~ 3 Gyr. The 'usual' 
complement of well-known old GCs (i.e., GCs of similar age 
as the majority of the Galactic GCs) is also present. 

As discussed in §2.6, to estimate the ages of our sample 
GCs accurately we required that our GC sample have both 
independently determine d metallicities and re d dening values. 
For m etalhc ity, we used 
(120001) . and iPerrett etal 



Huchra et alj ( Il991h . iBarmbv etal J 
(I2002h as our reference data set. 
Since all of these authors det ermined the M3 1 GC metal- 
licities using the calibration of iBrodie & Huchra] (1990), all 
three metallicity determinations are on the same [Fe/H] scale 
and there are no systematic offsets among any of these data 
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Fig. 4. — Best-fitting integrated SEDs of the galev SSP models shown in relation to the intrinsic SEDs for our sample GCs. The photometric data points are 
represented by the symbols with error bars (vertical error bars for uncertainties and horizontal ones for the approximate wavelength coverage of each filter). Open 
circles represent the calculated magnitude of the model SEDs for each filter. 
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Many M3 1 GCs are resolved in HST observa- 
tions. Some authors, includin g lAihar et alJ (Il996l). 



iFusi Pecci et"an (1 19961) 
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Fig. 5. — Age distribution of our sample GCs and GC candidates in M31. 



sets (see details in iPerrett et al.ll2002l) . However, individual 
metallicity differences exist, which may affect our age es- 
timates. Twelve GCs and GC candidates have two metal- 
licity determinations, of which B004, B219, and B238 ex- 
hibit the largest differences (> 0. 2 dex). The metal lici- 
ties of B004, B219, and B238 from IPerrett et al.1 (120021) are 
-0.31 + 0.74, -0.01 + 0.57, and -0.57 + 0.66 dex, com- 
pared to_2l;26± 0^9^0.53 + 0.53, and -1.22 + 0.76 dex 
from' Huchra et alJ (1199 Ih . Large metallicity differences lead 
to large age differences. The ages of B004, B219, and B238 
are estimated at 4.10 + 0.55, 2.50 + 0.15, and 5.00 + 0.45 
Gyr (based on the metalhcities of IPerrett et al.ll2002l) . com- 
pai-ed to 14.40 + 0.75, 11-60 + 1. 45 and 14.40 + 0.80 Gyr 
(metalhcities of Huc hra et al.lll99lh . This implies that the ac- 
curacy of the metallicity determinations is very important for 
the corresponding age estimates. The age differences for the 
other 10 GCs are less than 1 G yr based on th e metalhcities 
of both lPerrettetalJ (l2002h and lHuchraet all (1199 Ih . except 
forB045: 8 80+ 1 .45 Gyr versus 6.30 + 45 Gyr based on 
IPerrett et all (l2002h and lHuchra et alJ (ll991|) metalhcities. re- 
spectively. In general, the three different sources of spectro- 
scopic metalhcities provide homogeneous age estimates. For 
B004, B21 9, and B238 the sig nal- to-noise ratios of th e obser- 
vations of Huch ra et all (1 19911) and lPerrett et al.1 (l2002l) are too 
low. High-quality spectral observations of these three GCs are 
needed . In addition, we point ou t that the metallicity calibra- 
tion of iBrodie & Huchra) (1 19901) is solely based on old GCs. 
However, the sample GCs and GC candidates discussed in this 
paper are estimated to be young or of intermediate age, so that 
the age estimates may be somewhat biased by the adopted cal- 
ibration; 

iBarmbv et all (|2000|) discovered that M31 contains GCs 
exhibiting strong Balmer lines and A-type spectra, from 
which one infers that these GCs must be very young. 
iBeaslev et~a l. (2004) and Puziaetal. (200 5) confirme d this 
conclusion of Barmb v et all (120001) . iBurstein et all (12004) and 
iFusi PeccTet al. (2005|) in creased the sample of young M3 1 
GCs to 67. Very recently. ICaldwell et all (l2009l) determined 
the ages and reddening values of 140 young clusters in M31 
by comparing the observed spectra with model spectra, and 
these clusters are less than 2 Gyr old. Most have ages between 
10** and 10^ yr. The ages of the M31 clusters determined in 
this paper are in agreement with previous determinations. 



, , , .Richetal .' (1996"), 'Holland et aU 

Jablonka et all (120001) . IwTlhams & Hodge (200k 



[Wilhams & Hodgd (l2001b), and Rich et alJ (12005h . used 
WFPC2 images to construct CMDs for determination of the 
clusters' metalhcities, reddening values, and ages. However, 
these CMDs are usually not deep enough to show conspicuous 
MSTOs a nd thus be useful for robus t age determinations. In 
fact, on lvlWihiams & Ho d (1200 lab and lWihiam s & Hodgd 
dMl bl) managed to estimate the ages for four blue massive, 
compact star clusters and 79 candidate young star clusters by 
fitting isochrones to the stellar photometry. 

Our sample contains fou r GCs in common withlAihar et al] 
( 1996) (B006 and B045), Fusi Pecci et al. (1996") (also B006 
and B045), and Richetal. (2005) (B012 and B233). How- 
ever, only B012 (Rich et al. 2005) could be older than about 
8 Gyr (see Gallart et al. 2005), given the presence of a promi- 
nent blue horizontal branch, which compares rather poorly 
with the age obtained here (~ 2 Gyr). However, even with 
multi-passband photometry spanning from the FUV to Ks, we 
can only determine cluster ages in a statistical sense (also see 
iGaUagher & GrebellllOOllde Griis et al.ll2005h . 

This discrepancy of a few Gyr highlights the difficulty of 
obtaining age estimates of unresolved intermediate-age clus- 
ters based on multi-passband photometry, given that the color 
evolution of SSPs is only minimally age dependent once the 
population has reached an age of ~ 1-3 Gyr. In addi- 
tion, although the general age distribution of an entire clus- 
ter population (in a given galaxy) can be retrieved fairly self- 
consistently, the clusters' individual age determinations de- 
pe nd rather strongly o n the approach taken to fitting their ages 
(cf. lde Griis et alj200 5). This, combined with a strong depen- 
dence on the adopted reddening and metallicity, results in in- 
dividual age estimates of intermediate-age clusters associated 
with large uncertainties. 

Cluster ages can also be derived by compa ring observed 
with model spectra. Cr oss-i dentification o f iBeasl ev et alj 
(l2004 . lPuzia et al.l(l2005h . and lCaldwell etal. (2009) with the 
sample in this paper reveals that only six GCs (and GC can- 
didates) overlap with Puzia et al. (2005) (B006, B012, B045, 
and B232) and Caldwell et alj (2009) (B049 and B195). The 
age deerminati ons of two of these ( BO 12 and B232) are incon- 
sistent between [Puzia et aP (l2005l) a nd this paper: 10.2 + 2.9 
Gyr and 9.0 ± 3.3 Gyr obtained bv iPuzia et al] (l2005h com- 
pared to 2.0 ± 0. 1 Gyr and 2.0 ± 0. 1 Gyr, respectively, obtained 
in this paper. The ages of the other GCs are consistent among 
the different determinations. 

In fact, the ages of GCs derived by different authors based 
on a range of methods are not always consis tent. For ex- 
ample , the ages of B29 2 and B327 derived by iBeasley et alJ 
(2004) and 'Puzia et al.' ("2005) are 2.748 ± 1.151 Gyr and 
0.080 ± 0.929 Gyr CBeaslev et al...20 04) compared to 9.2 ±3.3 
Gyr and 5.4 ± 1 .4 Gyr (iPuzia et al.ll20 05). respectively. On the 
other hand. ICaldwell et alJ (1200 9*) estimated the age of B327 
at 0.050 Gyr. In addition, the ages of clusters derived in the 
same paper but based on different line-index measurements 
are not always consistent either and ma y indeed differ sig nifi- 
cantly (see the upper panels of Fig. 5 in lPuzia et alj|2005l) . 
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Table 1 

batc filter parameters and observational statistics 



FUter ^central (A) FWHM (A) N" Exp.'' rms'^ 



c 


4210 


320 


3 


01:00 


0.015 


d 


4540 


340 


3 


01:00 


0.009 


e 


4yzj 


ion 


D 


Ul.UU 


U.UID 


f 


5270 


340 


3 


01:00 


0.006 


g 


5795 


310 


3 


01:00 


0.003 


h 


6075 


310 


3 


01:00 


0.003 


i 


6656 


480 


3 


01:00 


0.003 


J 


7057 


300 


3 


01:00 


0.008 


k 


7546 


330 


3 


01:00 


0.004 


m 


8023 


260 


3 


01:00 


0.003 


n 


8480 


180 


6 


02:00 


0.004 





9182 


260 


6 


02:00 


0.003 


P 


9739 


270 


6 


02:00 


0.009 



" Number of exposures for each BATC filter 
* Exposure time (in hr:min) 
■^Zero-point error (in mag) 
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Table 2 

batc nvrrermediate-band photometry of our sample of 39 gcs and gc candidares in m31. 

Name c d e f g h i j k m n o p 

B004 17.71 17.49 17.22 17.06 16.72 16.61 16.44 16.34 16.21 16.03 16.09 15.91 16.01 

0.130 0.012 0.008 0.009 0.009 0.008 0.008 0.009 0.010 0.009 0.021 0.010 0.024 

B006 16.68 16.17 15.84 15.67 15.29 15.25 15.09 14.98 14.82 14.67 14.70 14.58 14.53 

0.009 0.005 0.004 0.004 0.004 0.004 0.004 0.003 0.004 0.003 0.008 0.005 0.008 

B008 17.83 17.33 17.04 16.88 16.51 16.43 16.26 16.14 16.01 15.87 15.73 15.80 15.73 

0.016 0.010 0.007 0.008 0.008 0.007 0.008 0.008 0.009 0.007 0.015 0.011 0.019 

BOlO 17.48 17.21 16.94 16.79 16.50 16.44 16.25 16.17 16.10 15.95 16.00 15.83 15.82 

0.014 0.010 0.007 0.008 0.008 0.007 0.008 0.008 0.011 0.008 0.019 0.012 0.023 

B012 15.90 15.58 15.34 15.20 14.91 14.85 14.71 14.63 14.55 14.40 14.45 14.36 14.35 

0.006 0.004 0.003 0.003 0.003 0.003 0.003 0.003 0.004 0.003 0.007 0.004 0.007 

BOB 18.35 17.86 17.51 17.36 16.99 16.92 16.75 16.64 16.47 16.33 16.36 16.24 16.19 

0.022 0.012 0.009 0.010 0.010 0.010 0.009 0.010 0.012 0.009 0.029 0.015 0.027 

B016 18.64 18.08 17.78 17.62 17.24 17.15 16.95 16.82 16.68 16.51 16.56 16.33 16.33 

0.027 0.016 0.011 0.012 0.012 0.011 0.011 0.010 0.014 0.010 0.028 0.016 0.029 

B019 15.75 15.45 15.15 14.99 14.61 14.53 14.38 14.25 14.11 13.96 13.89 13.80 13.76 

0.006 0.004 0.003 0.003 0.003 0.003 0.003 0.003 0.003 0.003 0.006 0.004 0.007 

B020 15.61 15.29 15.06 14.92 14.54 14.48 14.36 14.24 14.12 13.95 13.99 13.92 13.85 

0.006 0.003 0.003 0.003 0.003 0.003 0.003 0.003 0.003 0.003 0.006 0.004 0.007 

B022 17.96 17.74 17.53 17.40 17.09 17.05 16.96 16.89 16.80 16.65 16.65 16.59 16.61 

0.022 0.016 0.013 0.014 0.014 0.013 0.015 0.016 0.022 0.017 0.042 0.028 0.047 

B026 18.58 18.09 17.70 17.55 17.15 17.06 16.88 16.72 16.54 16.37 16.27 16.18 16.19 

0.031 0.023 0.016 0.018 0.020 0.015 0.017 0.015 0.018 0.018 0.033 0.022 0.035 

B035 18.45 17.98 17.68 17.55 17.17 17.06 16.92 16.81 16.62 16.49 16.52 16.39 16.33 

0.026 0.015 0.011 0.011 0.012 0.011 0.011 0.011 0.013 0.011 0.030 0.018 0.028 

B045 16.88 16.41 16.08 15.94 15.55 15.48 15.30 15.19 15.06 14.92 14.93 14.77 14.75 

0.010 0.006 0.004 0.005 0.005 0.004 0.004 0.004 0.005 0.004 0.009 0.007 0.010 

B047 18.36 17.95 17.68 17.59 17.23 17.15 17.02 16.96 16.86 16.69 16.77 16.70 16.73 

0.025 0.015 0.011 0.011 0.013 0.011 0.012 0.012 0.016 0.015 0.038 0.023 0.035 

B049 17.93 17.76 17.70 17.61 17.43 17.42 17.30 17.25 17.16 17.12 17.28 16.92 17.92 

0.025 0.022 0.020 0.021 0.024 0.019 0.021 0.024 0.026 0.021 0.061 0.040 0.412 

B050 17.68 17.32 17.03 16.87 16.49 16.42 16.27 16.14 16.01 15.84 15.78 15.70 15.73 

0.020 0.014 0.011 0.010 0.012 0.010 0.011 0.011 0.013 0.012 0.026 0.017 0.028 

B052 19.02 18.24 17.64 17.27 16.92 16.76 16.53 16.38 16.23 16.04 15.89 15.73 15.65 

0.046 0.020 0.012 0.011 0.012 0.010 0.009 0.009 0.012 0.009 0.020 0.013 0.021 

B062 18.86 18.17 17.66 17.24 16.96 16.74 16.52 16.44 16.31 16.13 15.98 15.83 15.79 

0.042 0.021 0.014 0.010 0.013 0.010 0.011 0.011 0.014 0.011 0.026 0.016 0.024 

B074 17.55 17.16 16.91 16.80 16.43 16.35 16.22 16.14 16.03 15.89 15.92 15.84 15.80 

0.014 0.010 0.007 0.007 0.008 0.007 0.008 0.008 0.009 0.008 0.027 0.013 0.020 

B081 17.05 17.01 16.87 16.81 16.59 16.42 16.44 16.18 16.08 16.01 15.95 15.53 15.86 

0.016 0.015 0.013 0.013 0.015 0.014 0.024 0.034 0.032 0.029 0.040 0.049 0.107 

B089 18.24 18.22 18.16 18.16 18.11 18.16 18.16 18.16 18.07 18.01 18.49 18.17 18.34 

0.027 0.026 0.026 0.024 0.039 0.034 0.041 0.045 0.050 0.047 0.126 0.090 0.163 
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Table 1 
Continued. 



Name c d e f 



BlOO 18.82 18.31 18.04 17.92 17.63 17.61 17.31 17.25 17.24 17.21 17.24 16.95 16.91 

0.050 0.030 0.026 0.026 0.027 0.025 0.026 0.048 0.031 0.033 0.065 0.037 0.170 

B129 18.41 17.71 16.83 16.59 16.07 15.79 15.47 15.21 14.99 14.71 14.62 

0.090 0.059 0.034 0.030 0.022 0.019 0.016 0.014 0.014 0.012 0.014 

B156 17.67 17.32 17.09 16.82 16.53 16.64 16.44 16.35 16.17 16.20 16.18 16.13 16.40 

0.018 0.011 0.008 0.008 0.008 0.010 0.012 0.014 0.012 0.031 0.021 0.035 0.008 

B168 19.50 18.76 18.23 17.96 17.29 17.13 16.78 16.60 16.28 16.07 15.95 15.80 15.69 

0.083 0.053 0.030 0.026 0.025 0.017 0.015 0.015 0.014 0.010 0.027 0.013 0.020 

B170 18.27 17.90 17.61 17.46 17.11 17.06 16.84 16.76 16.66 16.52 16.46 16.38 16.44 

0.025 0.020 0.012 0.013 0.014 0.012 0.012 0.014 0.015 0.014 0.033 0.018 0.141 

B195 18.97 18.78 18.61 18.57 18.38 18.36 18.35 18.11 18.20 18.04 17.96 17.96 ... 

0.046 0.035 0.028 0.035 0.044 0.033 0.044 0.044 0.062 0.059 0.109 0.074 ... 

B199 18.27 18.00 17.80 17.62 17.44 17.41 17.22 17.10 17.06 17.00 17.05 16.87 17.06 

0.024 0.018 0.012 0.013 0.016 0.013 0.014 0.016 0.018 0.017 0.039 0.022 0.060 

B207 18.04 17.74 17.53 17.36 17.14 17.10 16.93 16.84 16.81 16.73 16.71 16.59 16.64 

0.020 0.014 0.011 0.012 0.014 0.013 0.014 0.015 0.020 0.015 0.039 0.023 0.053 

B212 16.17 15.91 15.69 15.51 15.31 15.28 15.10 15.01 14.96 14.91 14.84 14.76 14.80 

0.007 0.005 0.004 0.004 0.004 0.004 0.004 0.004 0.005 0.004 0.011 0.006 0.012 

B219 17.25 16.90 16.63 16.44 16.15 16.05 15.88 15.76 15.62 15.46 15.46 15.32 15.33 

0.013 0.008 0.006 0.006 0.007 0.006 0.006 0.006 0.007 0.006 0.014 0.008 0.018 

B226 19.13 18.51 18.10 17.77 17.55 17.43 17.15 17.05 16.97 16.82 16.63 16.56 16.49 

0.044 0.029 0.015 0.015 0.018 0.015 0.016 0.016 0.022 0.016 0.041 0.021 0.045 

B230 16.44 16.33 16.13 15.98 15.75 15.68 15.56 15.47 15.41 15.21 15.21 15.20 15.20 

0.009 0.006 0.005 0.005 0.006 0.005 0.005 0.006 0.007 0.005 0.015 0.008 0.020 

B232 16.15 16.02 15.80 15.65 15.34 15.28 15.17 15.06 14.98 14.80 14.74 14.76 14.71 

0.008 0.005 0.004 0.005 0.005 0.005 0.005 0.005 0.006 0.006 0.011 0.008 0.016 

B233 16.41 16.19 15.94 15.81 15.44 15.38 15.26 15.15 15.01 14.85 14.85 14.73 14.90 

0.010 0.008 0.007 0.006 0.007 0.006 0.007 0.006 0.008 0.007 0.010 0.015 0.019 

B236 17.86 17.69 17.52 17.38 17.08 17.03 16.88 16.75 16.47 16.40 16.45 16.76 16.22 

0.021 0.017 0.012 0.012 0.015 0.012 0.012 0.018 0.023 0.040 0.021 0.039 0.164 

B237 17.92 17.70 17.45 17.31 17.01 16.94 16.79 16.70 16.63 16.48 16.50 16.42 16.47 

0.020 0.015 0.012 0.012 0.015 0.013 0.015 0.015 0.016 0.018 0.033 0.024 0.145 

B238 17.39 17.02 16.73 16.58 16.23 16.13 15.97 15.88 15.74 15.66 15.58 15.58 15.45 

0.014 0.007 0.006 0.006 0.007 0.012 0.007 0.016 0.016 0.014 0.019 0.011 0.055 

B239 17.88 17.79 17.51 17.43 17.01 16.90 16.77 16.66 16.52 16.43 16.33 16.32 16.15 

0.156 0.062 0.040 0.037 0.033 0.025 0.017 0.032 0.032 0.026 0.037 0.063 0.103 
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Table 3 

galex, broad-band, and 2mass photometry of the 39 m31 gcs and gc candmates. 



Name 


t 
c' 


njv 


NUV 


U 


B 


V 


R 


I 


J 


H 




B004 


1 




22.25 


18.29 


17.87 


16.95 


16.36 


15.73 


14.91 


14.36 


14.19 








0.07 


0.03 


0.01 


0.01 


0.02 


0.01 


0.02 


0.07 


0.05 


B006 


1 




21.41 


16.94 


16.49 


15.53 


14.97 


14.31 


13.48 


12.75 


12.61 








0.04 


0.02 


0.01 


0.01 


0.01 


0.01 


0.03 


0.03 


0.03 


BOOS 


1 




22.59 


18.16 


17.66 


16.56 


16.21 


15.51 


14.68 


14.17 


13.98 








0.12 


0.08 


0.05 


0.05 


0.05 


0.05 


0.05 


0.07 


0.06 


BOlO 


1 


21.93 


20.87 


17.65 


17.50 


16.66 


16.12 


15.48 


14.76 


14.41 


14.07 






0.08 


0.03 


0.02 


0.01 


0.01 


0.01 


0.01 


0.03 


0.07 


0.06 


B012 


1 


20.10 


19.02 


15.99 


15.86 


15.13 


14.62 


14.08 


13.36 


12.79 


12.72 






0.02 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0.03 


0.03 


0.03 


B013 


1 






18.56 


18.06 


17.19 


16.60 


15.96 


15.18 


14.61 


14.34 










0.05 


0.02 


0.01 


0.02 


0.02 


0.03 


0.07 


0.05 


B016 


1 






18.86 


18.58 


17.58 


16.85 


16.15 


15.18 


14.18 


14.05 










0.08 


0.04 


0.01 


0.03 


0.02 


0.08 


0.07 


0.10 


B019 


1 


20.81 


19.97 


16.36 


15.94 


14.93 


14.31 


13.74 


12.86 


12.10 


11.96 






0.04 


0.02 


0.01 


0.01 


0.01 


0.05 


0.01 


0.02 


0.03 


0.02 


B020 


1 


20.05 


19.20 


15.98 


15.74 


14.91 


14.37 


13.65 


12.97 


12.26 


12.21 






0.02 


0.01 


0.08 


0.05 


0.05 


0.05 


0.05 


0.02 


0.03 


0.03 


B022 


1 


22.69 


21.20 


18.14 


18.09 


17.36 


16.97 


16.35 


15.75 


15.04 


15.48 






0.16 


0.04 


0.08 


0.02 


0.01 


0.02 


0.02 


0.08 


0.10 


0.12 


B026 


1 






19.14 


18.60 


17.53 


16.88 


16.22 


15.10 


14.48 


14.00 










0.06 


0.02 


0.01 


0.03 


0.02 


0.08 


0.07 


0.05 


B035 


1 




22.61 


18.52 


18.37 


17.48 


16.81 


16.24 


15.30 


14.67 


14.42 








0.10 


0.08 


0.03 


0.01 


0.02 


0.02 


0.08 


0.07 


0.10 


B045 


1 




21.07 


17.09 


16.72 


15.78 


15.19 


14.54 


13.73 


13.00 


12.89 








0.03 


0.02 


0.01 


0.01 


0.01 


0.01 


0.03 


0.04 


0.03 


B047 


1 


22.68 


21.30 


18.32 


18.23 


17.51 


16.88 


16.30 


15.86 


15.24 


15.47 






0.15 


0.04 


0.06 


0.02 


0.01 


0.03 


0.02 


0.08 


0.10 


0.12 


B049 


1 




21.55 


18.26 


18.08 


17.56 


17.11 


16.87 


15.61 


15.33 


14.68 








0.06 


0.09 


0.04 


0.01 


0.04 


0.04 


0.08 


0.10 


0.10 


B050 


1 




22.18 


18.09 


17.76 


16.84 


16.27 


15.66 


14.72 


14.22 


13.96 








0.09 


0.05 


0.02 


0.01 


0.02 


0.01 


0.03 


0.07 


0.05 


B052 


4 






19.80 


18.62 


17.21 


16.54 


15.77 


14.70 


14.01 


13.40 










0.08 


0.02 


0.01 


0.02 


0.02 


0.04 


0.07 


0.05 


B062 


4 






19.33 


18.58 


17.24 


16.61 


15.82 


14.89 


14.21 


13.66 










0.08 


0.02 


0.01 


0.02 


0.01 


0.04 


0.07 


0.05 


B074 


1 


22.12 20.75 


17.54 


17.40 


16.65 


16.14 


15.58 


14.83 


13.95 


14.11 






0.07 


0.02 


0.03 


0.01 


0.01 


0.01 


0.01 


0.02 


0.04 


0.04 


B081 


1 


21.73 


20.47 


17.60 


17.34 


16.80 


16.36 


15.73 


14.82 


14.01 


13.96 






0.13 


0.04 


0.02 


0.01 


0.01 


0.02 


0.02 


0.03 


0.07 


0.05 


B089 


2 


19.89 


19.62 


17.96 


18.28 


18.18 


18.22 


17.70 












0.03 


0.02 


0.05 


0.04 


0.03 


0.06 


0.06 
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Table 3 
Continued. 



Name 


c 


FUV 


NUV 


U 


B 


V 


R 


/ 


J 


H 




BlOO 


1 


22.37 




18.94 


19.05 


17.91 




17.77 


15.85 


14.70 


14.67 






0.19 




0.08 


0.07 


0.05 




0.07 


0.08 


0.07 


0.10 


B129 


1 








19.56 


17.40 




14.69 


13.25 


12.40 


12.19 












0.05 


0.05 




0.05 


0.03 


0.03 


0.04 


B156 


1 




20.99 


17.89 


17.63 


16.84 


16.37 


15.87 














0.09 


0.02 


0.02 


0.01 


0.02 


0.05 








B168 


1 






20.82 


19.23 


17.63 


16.69 


15.72 


14.52 


13.43 


13.37 










0.08 


0.06 


0.01 


0.03 


0.02 


0.06 


0.04 


0.08 


B170 


1 






18.90 


18.37 


17.39 


16.80 


16.17 


15.38 


14.75 


14.61 










0.06 


0.02 


0.01 


0.02 


0.02 


0.08 


0.07 


0.10 


B195 


2 






19.94 


18.97 


18.57 


18.02 


17.59 
















0.08 


0.06 


0.01 


0.07 


0.05 








B199 


1 




21.53 


18.45 


18.37 


17.60 


17.03 


16.57 


16.06 


15.42 


15.39 








0.10 


0.08 


0.03 


0.01 


0.03 


0.02 


0.10 


0.10 


0.12 


B207 


1 


21.64 


21.04 


18.26 


18.07 


17.33 


16.81 


16.33 


15.67 


14.42 


14.78 






0.11 


0.06 


0.03 


0.02 


0.01 


0.02 


0.02 


0.05 


0.07 


0.08 


B212 


1 


20.27 


19.16 


16.23 


16.22 


15.48 


15.00 


14.48 


13.82 


13.17 


13.11 






0.05 


0.02 


0.02 


0.01 


0.01 


0.01 


0.01 


0.03 


0.04 


0.05 


B219 


1 




21.59 


17.74 


17.32 


16.39 


15.82 


15.19 


14.32 


13.71 


13.51 








0.10 


0.08 


0.05 


0.05 


0.05 


0.05 


0.02 


0.04 


0.04 


B226 


2 22.09 


21.61 


19.08 


19.04 


17.65 




16.32 


15.21 


14.47 


14.14 






0.14 


0.09 


0.08 


0.05 


0.05 




0.05 


0.08 


0.07 


0.10 


B230 


1 


20.69 


19.46 


16.78 


16.77 


16.05 


15.61 


15.13 


14.43 


13.92 


13.85 






0.08 


0.03 


0.02 


0.01 


0.01 


0.01 


0.01 


0.02 


0.04 


0.05 


B232 


1 


20.67 


19.49 


16.53 


16.38 


15.70 


15.20 


14.65 


13.94 


13.36 


13.25 






0.06 


0.02 


0.01 


0.01 


0.01 


0.01 


0.01 


0.03 


0.04 


0.05 


B233 


1 


21.27 


20.04 


16.82 


16.61 


15.80 


15.27 


14.76 


13.90 


13.32 


13.21 






0.05 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0.02 


0.04 


0.03 


B236 


1 




21.07 


18.21 


18.20 


17.38 


16.97 


16.24 














0.07 


0.08 


0.02 


0.01 


0.02 


0.02 








B237 


1 




21.25 


18.03 


17.87 


17.10 


16.57 


16.05 


15.47 


15.06 


14.91 








0.05 


0.02 


0.02 


0.01 


0.02 


0.02 


0.04 


0.10 


0.09 


B238 


1 




21.91 


17.73 


17.39 


16.42 


15.86 


15.22 


14.46 


13.72 


13.67 








0.08 


0.02 


0.01 


0.01 


0.01 


0.01 


0.04 


0.04 


0.05 


B239 


1 






18.49 


18.10 


17.08 


16.65 


16.09 


15.31 


14.55 


14.55 










0.04 


0.02 


0.01 


0.02 


0.02 


0.04 


0.07 


0.08 



t New classification flag, following RBC V3.5 notation. 1 = confirmed GC, 2 = GC candidate, 4 = confirmed galaxy 
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Table 4 

Reddening values and metallicities for our 39 M31 GCs and GC candidates. 



Name 


£L{n — V ) 


lei. 


[te/ti] 


ret. 


B004 


U.U/± 0.02 


1 


A T 1 r n 1 1 

— 0.31± 0.74 


1 


B006 


U.Uy± (J.Uz 


1 


A C O _i_ A 1/1 


1 


BOOS 


0.21 


2 


— 0.41± 0.38 


1 


BO 10 


0.2i± 0.01 


1 


— 1.77± 0.14 


1 


B012 


0. 12± 0.01 


1 


-1.65± 0.19 


1 


B013 


0.1j± 0.02 


1 


1 Ai 1 r\ An 

— 1.01± 0.49 


1 


B016 


A OA 1 A A1 

U.3U± 0.02 


1 


A "70 1 A 1 A 

—0. /8± 0.19 


1 


B019 


A lA 1 A A 1 

0.20± 0.01 


1 


1 AA 1 A AT 

— 1.09± 0.02 


1 


BUZU 


n 1 A A 1 


1 


1 r\n _i_ A 1 A 
— i.(J/± U.IU 


3 


B022 


A A /I 1 A Af 

0.04± 0.03 


1 


1 ^ A 1 A AT 

— l.o4± U.U/ 


4 


B026 


A 1 C 1 A A1 

0.15± 0.02 


1 


/I A 1 1 A T O 

0.01± 0.38 


1 


B035 


A T7 1 A AC 

0.2 /± 0.05 


2 


A TA 1 A C /I 

— 0.20± 0.54 


1 


B045 


A 1 O 1 A /"\ 1 

0.i8± 0.01 


1 


— 1.05± 0.25 


1 


B047 


u.uy± U.UZ 


1 


— i.OZ± U.4i 


1 


B049 


r\ t /I 1 A Ai 

0.1o± 0.02 


1 


— 2.14± 0.55 


1 


B050 


A 1 ^ 1 A A1 

0.24± 0.01 


1 


1 /I T 1 A IT 

— 1.42± 0.3 / 


1 


B052 


A TI 1 A f\A 

0.23± 0.04 


1 


/I 1 O 1 AIT 

U.12± 0.1 / 


4 




U.zo± U.OJ 


1 


A /I T-L A 1 1 
— U.4/± U. i 1 


4 


B074 


A 1 A 1 A A 1 

o.iy± 0.01 


1 


1 O O 1 A Aj^ 

-1.88± U.Uo 


1 


B081 


A 1 1 _i_ A AT 

U. i i± U.Uz 


1 


1 "7/1 _i_ A /I A 

— i. /4± U.4U 


1 


B089 










BlOO 


0.48± 0.08 


1 


O O 1 A 1 A 

-2.21± 0.10 


4 


B129 


1.1 6± 0.06 


1 


1 1 1 1 ATT 

— 1.21± U.32 


1 


B156 


A "1 A 1 A A1 

0.1 0± 0.02 


1 


-1.51± 0.38 


1 


1) 1 y-o 

Bloo 


A Cyl _i_ A AC 

0.54± O.OD 


1 


A 1 ATI 
— U.iZ± U.Zi 


4 


B170 


A 1 A_i_ A AT 

0.10± 0.02 


1 


A C /I _i_ A O /I 

— U.54± U.z4 


1 


B195 


Alii A AA 

0.12± 0.00 


1 


1 AO 1 A 

— 1.48± U.o3 


4 


B199 


A 1 A_i_ A AO 
0.10± 0.02 


1 


— 1.59± 0. 1 1 


1 


B207 


A AC 1 A AT 

0.05± 0.02 


2 


A O 1 1 A CA 

— 0.81± 0.59 


1 




U. i j± U.Ui 


1 
i 


1 7S-t- n 1 

i . / J± U. 1 J 


J 


B219 


0.05± 0.03 


1 


-0.01± 0.57 


1 


B226 


1.08± 0.06 


1 






B230 


0.15±0.01 


1 


-2.17± 0.16 


1 


B232 


0.14± 0.01 


1 


-1.83± 0.14 


1 


B233 


0.17± 0.01 


1 


-1.59± 0.32 


3 


B236 


0.07± 0.05 


1 


-1.01± 0.17 


4 


B237 


0.14± 0.02 


1 


-2.09± 0.28 


1 


B238 


0.11±0.02 


1 


-0.57± 0.66 


1 


B239 


0.09± 0.01 


1 


-1.18± 0.61 


2 



"The reddening values are taken from'Fan et al.' ("2008^ (ref.=l) and'Barmbv et al.' FxM) (ref.=2). 
''The metallicities are taken from .Perrett et aH2Q0Z^ (ref.=l). ,Bai-mbv et al. t200a) (ref.=2). iHuchia et al. (1991.) (ref.=3), and lFan et all j2008h (ref.=4). 
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Table 5 

Ages estimates for 35 GCs and GC candidates in M31. 



Name 


Age 




Name 


Age 






(uyr) 


(per degree of freedom) 




(Cryr) 


(per degree of freedom) 


dUU4 


A 1 _i_ fl <C 

4.1U ± \J.DJ 


J. 13 


D 1 An 
BlOO 


A CA _i_ A 1 A 
O.jO ± O.IU 


1 A 1Q 
14.J0 


dUUo 


12. DU ± U.OJ 


1 in 


D 1 on 


1 C 1 A _i_ A TA 
ID. 10 ± O./O 


9.0U 


BOOS 


2. 00 ± 0.10 


6.54 


B156 


A AA 1 A 

4.90 ± 0.65 


2.09 


BOlO 


1 OA _i_ A 1 A 
l.oU ± U.IK) 


1 .08 


B168 


12. oO ± O.zO 


3.27 


B012 


r\f\ 1 A 1 A 

2.00 ± 0.10 


1.54 


B170 


4.00 ± 0.45 


1.39 


B013 


12.00 ± 2.00 


0.96 


B195 


CI "7A 1 ("lie 

0. /O ± 0.15 


0.82 


Dm ^ 
BUlo 


1 Ar\ _i_ A ^A 

2.40 ± 0.30 


1 

i.y5 


B199 


Q 1A 1 A CC 

j.jO ± 0.55 


1.38 


B019 


O 1 A 1 A 1 A 
2.10 ± 0.10 


12.01 


B207 


1 OA 1 A 1 A 
1.20 ± 0.10 


7.99 


B020 


1 Rft -1- n 1 ft 




B212 


1 sn -1- ft 1 n 


ft 7S 


B022 


3.40 ±0.15 


4.11 


B219 


2.50 ±0.15 


3.06 


B026 


3.50 ± 0.25 


3.28 


B230 


1.60 ±0.10 


4.59 


B035 


1.00 ±0.10 


3.97 


B232 


2.00 ±0.10 


2.62 


B045 


8.80 ± 1.45 


0.78 


B233 


2.30 ±0.10 


3.73 


B047 


2.80 ± 0.20 


4.22 


B236 


2.00 ± 0.25 


2.63 


B049 


1.60 ±0.10 


7.82 


B237 


3.50 ± 0.35 


1.41 


B050 


16.00 ±0.30 


2.12 


B238 


5.00 ± 0.45 


2.01 


B074 


2.10 ±0.15 


2.12 


B239 


14.50 ± 2.05 


1.70 


B081 


2.10 ±0.20 


7.82 









